Communications Physics

Maintaining a prescribed composition in engineered microbial consortia is difficult because small fitness differences can drive competitive exclusion. We study a two-strain consortium in continuous culture and develop a feedback architecture that regulates composition by selectively slowing the fast strain as a function of the population ratio. At the population level, we derive an idealized ratio-feedback law with a tunable positive coexistence equilibrium. We then propose a biomolecular realization using orthogonal quorum sensing, an sRNA-based ratiometric controller, and a ppGpp-mediated growth actuator. Exploiting the separation between slow population growth and faster intracellular controller dynamics, we use singular perturbation theory to show that, for sufficiently fast controller dynamics, the full implementation model inherits the coexistence equilibrium and its local stability properties from the reduced model. Numerical simulations validate the reduction and show how weaker timescale separation or loss of the assumed molecular regime degrades performance.

Radiation-induced DNA double-strand breaks (DSBs) drive cellular mortality, mutagenesis, and severe evolutionary bottlenecks. While classical phenomenological models, such as the Linear-Quadratic (LQ) framework, reliably predict macroscopic population survival, they obscure the intrinsic single-cell stochasticity that governs critical rare events like tumor recurrence or the emergence of radioresistant persisters. To bridge this divide, we develop a mathematically exact stochastic differential equation (SDE) framework that models continuous DSB induction and repair as a Feller square-root process. By deriving exact closed-form expressions for the foci moments, we establish a highly efficient Maximum Likelihood Estimation (MLE) pipeline that circumvents computationally exhaustive Monte Carlo simulations, allowing the direct extraction of deterministic repair velocities and intrinsic molecular noise from empirical single-cell $\gamma$-H2AX data. Integrating this kinetic model with a cumulative damage hazard via the Feynman-Kac formalism, our framework seamlessly recovers the classic macroscopic LQ survival topology from microscopic first principles. Furthermore, systematic sensitivity analysis uncovers a fundamental evolutionary duality: while initial physical damage operates additively, ultimate cellular fate is driven by a nonlinear survival response governed by the trade-off between the damage hazard rate and intrinsic molecular noise strength. Crucially, we demonstrate that this molecular noise inherently enhances population survival. Governed by Jensen's inequality, stochastic variance acts as a non-genetic bet-hedging mechanism that buffers the population by favoring cells with transiently low damage loads. Ultimately, this exact stochastic framework bridges microscopic biophysics and macroscopic demographics, offering deep mechanistic insights into the evolutionary roots of radioresistance.

Many biological processes rely on mechanical forces, with protein molecules acting as key mediators. Understanding how proteins respond to mechanical stress is essential for conditions including cardiomyopathy and muscular dystrophy. Natural proteins such as dystrophin and utrophin are composed of heterogeneous folding domains with distinct mechanical properties; deciphering domain-level behavior provides insights into disease mechanisms and informs therapeutic strategies. Single-molecule force spectroscopy (SMFS) enables probing the mechanical properties of entire proteins, yet current approaches struggle to identify heterogeneous folding domains, particularly without prior knowledge. Here, we present the first automated framework to identify heterogeneous folding domains in SMFS data, applying both existing clustering methods and a novel physics-aware deep clustering architecture, LatentUnfold. LatentUnfold learns complementary latent representations from force magnitude and the force-extension physical relationship through dual autoencoders, jointly optimized for clustering assignments. We apply our framework to experimental SMFS data collected from a synthetic two-domain protein (ddFLN4-Titin I27) as well as natural protein constructs of dystrophin and utrophin, with Monte Carlo simulated datasets serving as controlled validation. For the synthetic protein, we recover mechanical properties consistent with previously reported values for each domain. For the natural proteins, we uncover two mechanically distinct domain populations - corresponding to the N-terminal domain and spectrin-like repeats - with differences in both unfolding force and contour length increase, and reveal different unfolding order between them for the first time. This work enables domain-level biological inference, overcoming prior limitations that relied on averaging and overlooked heterogeneity, thus advancing the understanding of mechanical behavior in protein unfolding.

Biological aging is accompanied by systematic changes in epigenetic modifications and chromatin organization. The reversal of the effects of aging, rejuvenation, is experimentally achieved by the transient induction of factors that modify these marks in cells and organisms. Here, we show that key features of rejuvenation experiments emerge from the biophysical interplay between dynamic epigenetic marks and the three-dimensional conformation of chromatin. Using a minimal field theory and molecular dynamics simulations, we show that the system responds in three distinct temporal regimes. The intermediary regime fulfills necessary conditions for successful rejuvenation. In this regime, the system spends time near a separatrix, allowing for high epigenetic plasticity, while memory retained in the chromatin conformation enables restoration of the original epigenetic correlations. Analysis of sequencing data further supports the predicted coupling between chromatin compaction and epigenetic correlations. Our results provide a physical explanation for how rejuvenation may remodel age-associated epigenetic states without irreversibly erasing cellular identity. We identify a general mechanism by which memory stored in a slow structural variable permits reversible remodeling of a faster internal state.

In linear polysomes, excluded-volume interactions among ribosomes can induce dimensional reduction of mRNA. Yet linear architectures allow steric stress to relax at open ends-- limiting how strongly crowding can remodel the mRNA's structure and dynamics. Using coarse-grained molecular-dynamics simulations, we compare circular and linear polysomes over a range of ribosome densities. Circular closure selects a predominantly quasi-planar global conformational ensemble, as indicated by a shape dimensionality dshape {approx} 2 over a range of ribosome densities. Crucially, circular topology and ribosome crowding act cooperatively to suppress structural fluctuations. While closure alone or linear crowding reduces relative global size fluctuations ({Delta}Rg/Rg) only to {approx} 0.16, their combined effect drives this fluctuation down to {approx} 0.07. Within this stabilized architecture, increasing ribosome density drives a distinct in-plane reorganization: the ring becomes more isotropic, global size fluctuations are strongly suppressed, and the scaling exponent increases toward {nu} [~=] 0.74 - 0.77, consistent with two-dimensional self-avoiding walk-like value over the accessible finite-size window, 1000 [≤] N [≤] 4969. Closure shortens the radius-of-gyration decorrelation time of circular polysomes by 40-fold relative to matched linear systems, reflecting the topological elimination of free ends. Within this closureselected ensemble, ribosome crowding further reduces the decorrelation time by up to 20% at the highest density. A fluctuation-informed crossover model links the density dependence of the global scaling exponent to inter-ribosomal subchain statistics. These results distinguish the geometric role of circular closure from the density-dependent steric response that it enables, revealing a confined yet dynamically responsive conformational regime for circular polysomes.

The dynamics of gene regulatory networks are governed by intrinsic noise, stemming from the random nature of biochemical reactions, and by extrinsic noise, arising from fluctuations in cellular components and environmental conditions. Together, these sources can compromise the reliability of predictive computational models if not properly accounted for, and capturing both effects within a single framework remains a non-trivial task in computational biology. In this work, we propose an uncertainty quantification framework that addresses these two contributions jointly: intrinsic stochasticity is described through a partial integro-differential equation (PIDE) for the protein probability density function, whereas extrinsic noise is represented as parametric uncertainty in the kinetic parameters. The propagation of the uncertainty is carried out via an intrusive polynomial chaos expansion (PCE), in which the PCE coefficients are obtained from a stochastic Galerkin projection of the PIDE, yielding a coupled deterministic system that is solved with standard numerical methods. We illustrate the approach on a positive autoregulatory gene network with one and two uncertain kinetic parameters. The proposed approach accurately reproduces the mean, variance, and full protein probability density function, including the bimodal distributions, at a substantially lower computational cost.

Focal adhesions maintain force-bearing attachment between cells and the extracellular matrix but can also undergo dynamic remodeling. Their assembly and actomyosin tension are coupled through mechanochemical feedback. The processes underlying this feedback are not instantaneous and therefore involve a time delay. However, how this delayed feedback gives rise to stable adhesion maintenance or dynamic remodeling remains unclear. Here, paired time-lapse measurements of vinculin fluorescence and traction stress revealed distinct local adhesion-force dynamics, including low-fluctuation and recurrent fluctuation patterns. To examine how these patterns could arise, we formulated a minimal mechanochemical model coupling focal adhesion assembly and actomyosin force through delayed reciprocal feedback. The model exhibited stable and oscillatory modes depending on feedback strength, the balance of opposing feedback effects, and the effective feedback delay. Bistability and hysteretic switching also occurred in a subset of parameter space, and the oscillation period followed a power-law relation with the delay. These results suggest that stable adhesion maintenance and dynamic remodeling can emerge from a common mechanochemical feedback architecture.

Adhesion of cells is the key factor determining functioning of multicellular organisms. Viscoelastic properties of cells can be studied by multiple methods. However, attractiveness of cells or extracellular matrix without the elastic component (dispersive adhesion) is not accessible. We present an extension of force spectrometry technology: the Multivalent Adhesive Probe Atomic Force Microscopy (MAPA) that delivers dispersive adhesion maps of live cells and biosurfaces, and identifies differences unresolved by viscoelastic probing.

Modern medicine implicitly assumes that physiological responses to intervention are predictably determined by administered treatments. However, physiological systems containing intrinsic delays between the detection of a stimulus and the biological response may violate this assumption. We investigate the human glucose-insulin system as described by the Ultradian model and mathematically demonstrate that clinically relevant forcing protocols-such as pulsatile insulin delivery and step-wise glucose infusion, both commonly used in intensive care units (ICUs)-can induce sustained temporal chaos that may hamper accurate prediction of the physiological response. If not accounted for, these chaotic dynamics could create difficulties in achieving optimal dosing and timing when administering glucose and insulin in clinical or home care settings. This phenomenon, termed delay-induced uncertainty (DIU), arises from the interaction between physiological delay, intrinsic shear near a limit cycle, and external forcing. Using the Ultradian glucose-insulin model, we compute top Lyapunov exponents to quantify predictability. Across a range of pulsatile and step-wise forcing regimes, including stochastic amplitudes drawn from Markov processes, we observe positive Lyapunov exponents, indicating sustained chaos. Our results suggest that delayed endocrine regulation may fundamentally limit the predictive value of the models used to develop glycemic management strategies, with implications for clinical protocols in the ICU.

SynNotch-CAR circuits improve T cell specificity by coupling antigen recognition to inducible CAR expression. However, basal CAR expression without receptor activation, termed here as leakiness, can reduce the separation between killing of intended target cells and sparing of antigen-positive off-target cells, limiting target-cell discrimination. Here, we systematically quantified basal CAR expression for several synNotch-CAR designs and developed a coupled ordinary differential equation model to show that discrimination depends on basal output, CAR potency, and effector-to-target ratio. We introduced C-terminal tags such as fluorescent proteins, degron domains, endocytosis signals, and endoplasmic reticulum retention motifs as a strategy to reduce CAR leakiness. We found that fluorescent proteins and degron-containing tags reduced basal CAR surface expression while preserving antigen-induced CAR expression, improving discrimination of antigen-density sensing and combinatorial circuits in vitro. In xenograft models, fluorescent protein-tagged CARs improved discrimination by reducing activity against off-target cells while retaining activity against high-antigen tumors. Degron-containing constructs reduced basal CAR expression in vitro but showed suboptimal performance in vivo, revealing a trade-off between basal CAR suppression and induced CAR persistence. Together, these findings demonstrate that basal output expression is a key parameter for inducible genetic circuit designs and establish layered transcriptional and post-translational regulation as a strategy to improve the fidelity of inducible T cell circuits.

Spatial-omics workflows enable molecular analysis within tissue spatial context. Despite the prognostic value of tissue stiffness, these approaches have not incorporated direct, mechanical measurements. This omission reflects several challenges, including sample requirements, low throughput, specialized equipment, and complex data registration. Here, we introduce mechanoMaST (mechanics mapped to spatial transcriptomics), the first workflow to combine absolute mechanical measurements with spatial-omics. It pairs atomic force microscopy-based nanoindentation stiffness maps with spatial transcriptomics maps from adjacent tissue cryosections. The two modalities are then computationally co-registered to enable direct spatial correlation at 100 um resolution, with mapping accuracy quantified through error propagation, providing ground-truth mechanical data directly linked to spatial gene expression. We demonstrate mechanoMaST in human colorectal cancer liver metastasis, generating a spatial resource from 10 patients and revealing a four-gene stiffness signature. mechanoMaST is readily adaptable to other tissues across development and disease, and extendable to additional spatial-omics modalities in adjacent sections.

The shape of density-dependence governs species persistence, and ecosystem stability. Yet, whether per-capita growth declines sublinearily, or superlinearily with density remains hotly debated. Growth rates across the tree of life have been shown to decline sublinearly with density, whereas theory founded on resource competition predicts the opposite. Here, we resolve this discrepancy and show that sublinearity can readily emerge from geometric constraints on consumer interactions. By linking inter individual spacing, movement and interference rates, we derive two limiting-interference regimes, one of which the well-mixed limit recovers the form of classic Beddington DeAngelis interference response. We then developed an individual-based model from first principles which reproduces the derived sublinearity response, and further use empirical data from published consumer-resource experiments that also bears the signature of sublinear density-dependence. Further, embedding the interference mechanisms underlying the emergence of sublinear density-dependence in coexistence theory opens a new regime for species coexistence where classical theory fails to predict. Our framework indicates that non-consumptive interactions are not merely a correction to resource competition but might be a distinct axis along which diverse communities may potentially coexist.

Cell motility drives spatial pattern formation across diverse biological systems. Here, we engineer Escherichia coli motility in semi-solid agar to control Voronoi patterns in two and three dimensions, partitioning space into regions closest to their respective inoculation seeds. Consistent with our reaction-diffusion model, we observed that collisions between expansion fronts generate either biomass depletion (''gaps'') or accumulation (''anti-gaps''), governed by the relative diffusion rates of bacteria and nutrients. By engineering strains with distinct expansion rates and tuneable motility, and by integrating these experimental data into a dynamic Voronoi model, we achieved precise control over pattern geometry. This enabled the generation of gaps with varying widths, curved boundaries, asymmetric structures, seedless regions, and complex composite patterns. Together, these findings establish bacterial Voronoi patterns as a programmable platform for engineering multicellular spatial organization, with potential applications in synthetic biology and materials science.

Label-free Raman imaging of living cells is photon limited: at exposures compatible with cellular dynamics, single-pixel spectra carry about one count per channel on a dominant smooth background. We present an unmixing framework in which the decoder of a physics-constrained autoencoder is restricted to a data-driven spectroscopic dictionary: band centers,widths, and pseudo-Voigt shapes are measured from the dataset and fixed, and the network learns only nonnegative band amplitudes, a smooth B-spline background, and a per-pixel gain.First, on slit-scanning images of HeLa cells (532 nm) the dictionary yields spike-free component spectra that read as band tables, including a resonance-enhanced cytochrome-c-associated component matching literature spectra, and the most stable decomposition against the component number. Second, the dictionary and initialization calibrated at 1 s exposure perline transfer to 100 ms per line (12 s sweeps): cytochrome-c spectral identity survives a single sweep (correlation 0.92) while its map remains photon limited; the dictionary provides spectral physicality, and the transferred initialization prevents a structural collapse that global map correlations miss; in a measurement-derived phantom the dictionary estimator holds thecytochrome-c spectrum to 17-19{degrees} spectral angle at 100 ms, where classical factorizations and free decoders lose it (55-64{degrees}). Estimation on the count-equivalent detector output uses a calibrated shifted-Poisson quasi-likelihood. Third, evaluation must be time matched:correlation against a separately acquired reference saturates through slow specimen drift and acquisition mismatch rather than photon noise, and the self-consistency of learned denoisers is inflated by shared bias; time-matched self-consistency and independent cross-checks areproposed.

Delamination within stratified epithelia like the skin epidermis describes the detachment and upward motion of cells originating from the basal layer. Despite its fundamental importance for tissue development, homeostatic regeneration and repair, the mechanisms that drive delamination remain a longstanding open question. Upward motion follows cell shape changes, which are inherently driven by physical forces, but their role is elusive. Here, we investigate delamination in stratifying keratinocytes by combining imaging, force measurements and theoretical modeling. We identify a local change in force balance between differentiating cells and their environment as the key step initiating delamination. Within a homogeneous cell layer with apically polarized contractility, differentiation leads to actomyosin remodeling, redistributing cellular force exertion to the basal side. Such mechanical heterogeneity then results in differentiating cells experiencing and inward basal and outward apical forces that manifest in the formation of a +1 force defect and promote shape changes culminating in upward motion. Simultaneously, delaminating cells actively pull on their underlying neighbors, generating convergent tissue flows which close the basal layer below. Together, we propose a general physical description of delamination initiation, which may act across various multilayered epithelia.

A dental AI system and a dentist reading the same radiographs form a paired comparison. Published comparative studies often report the two arms separately against a reference standard, leaving the joint pattern of correctness between them unavailable for secondary paired inference. We show what that omission costs. The accuracy difference remains exactly identified; its sampling variance does not, so the report contains the estimate and not its uncertainty. On a study of 282 units, two published accuracies are consistent with 38 distinct joint tables whose confidence intervals differ in width by a factor of 2.5. The consequence is a three-zone decision map rather than a single threshold: differences at or below 1.06 points are non-significant under every compatible table, differences at or above 6.03 points are significant under every compatible table, and in between the published numbers cannot decide. We then show the omission is repairable at negligible cost. One additional integer, the number of units both arms classify correctly, identifies the joint table exactly and restores standard paired inference. For a panel of readers the pairwise dependences must arise from one joint distribution, a constraint that binds once three readers are present; publishing each reader's joint-correct count against a single reference reader cannot widen and may tighten every pairwise bound, and in a 7-arm experiment reduced them by a median of 37% even for pairs excluding that reference. Where the integer was never published we give DentalPair-Cert, an interval with finite-sample coverage uniformly over every admissible within-unit AI-dentist dependence under the independent-sampling-unit model, certified in both the nuisance maximization and the inversion. Across 4,200,000 simulated comparisons an independence analysis falls to 74.5% coverage with 12.2% type-I error; in a purposive sample of 9 recent comparative studies, 1 reported a paired test on discordant units.

Pouring a drink feels deceptively trivial, yet it requires guiding a boundary-free fluid into a vessel without spilling, overflowing, or toppling it -- a task at which robots remain notoriously brittle. How humans achieve this so effortlessly is unknown, as motor control has predominantly been studied in brief, highly constrained laboratory tasks, leaving the control principles underlying ecological tasks largely unknown. Here we measured continuous sensorimotor control during liquid pouring across various containers, vessels, and speed demands. Despite substantial variation in movement trajectories and durations, individuals maintained a strikingly invariant preferred fill level. Counterintuitively, fill level variability decreased at higher fill levels, and precision was maintained even under time pressure. A stochastic optimal control model combining a data-driven nonlinear approximation of flow dynamics with a cost that balanced individualised fill level, energy expenditure and flow-rate reproduced the behaviour. Humans thus pour optimally, given their sensorimotor limits and idiosyncratic notion of "full".

Thermodynamic models quantitatively describe interactions between transcription machinery and bacterial promoters. Contrary to conventional understanding, model analysis by Parisutham et al. (2025) attributes transcriptional inhibition by repressors to overstabilization of the RNA polymerase-promoter complex rather than prevention of its formation. Moreover, it suggests an inverse scaling relationship between basal promoter strength and transcriptional fold change, applicable to both repressor- and activator-mediated regulation. To reevaluate findings from this study, we systematically analyze empirical data and compare its framework with conventional thermodynamic models. In contrast to the inverse scaling relationship, data across multiple sources exhibit a peaked tradeoff between basal promoter strength and fold change, underscoring the importance of broad data coverage in revealing the full pattern required for reliable model inference. Furthermore, we identify the model assumption responsible for the apparent inverse scaling and misinterpretation of regulatory mechanisms. Relaxing this assumption enables the model to capture the peaked tradeoff and yield inferences consistent with established mechanisms of transcriptional repression and activation. We further derive a mathematical solution that connects basal expression to fold change for both repressor- and activator-regulated promoters. Our results underscore the importance of broad data coverage to avoid a blind-men-and-elephant interpretation and establish basal promoter strength as a key design parameter governing transcriptional regulation.

The folding thermodynamics of proteins are dominated by two opposing forces, the loss in backbone entropy and the packing of hydrophobic groups. The same forces are major contributors to the extension thermodynamics of elastic proteins with the distinction that both processes act in concert, favoring the higher chain and solvent entropy of a relaxed conformation. The relative entropic contributions specify the recoil mechanism; human elastin recoil is primarily driven by hydrophobic forces, whereas fly resilin has a rubber-like mechanism driven by backbone entropy. Despite the importance of elastic proteins to tissue biomechanics, few have been identified, let alone characterized to the same extent as elastin and resilin. We develop a thermodynamic framework that maps proteins by sequence-derived estimates of extension-induced backbone and solvent entropy changes. Putative elastic proteins are proposed and classified by recoil mechanism based on estimated thermodynamic features. Proteins that map to elastic regions are overrepresented by the skin proteome. The set of predicted elastic domains is further extended by incorporating sequence context embedded in protein language models. Protein domains with distinct thermodynamic recoil mechanisms cluster on the latent space manifold. Some of these domains are anticipated to have roles within molecular machines, expanding the scope of elastic protein function beyond mechanical materials like elastin and resilin.

Protein kinases function as dynamic, mechanically coupled nodes, yet the conformational drivers of multimeric activation remain unclear. Here, we present AlloQuant, a computational suite that translates AlphaFold3 structural ensembles into quantitative metrics of kinase regulation, including internal network rigidity, metastable-state populations, and sub-angstrom conformational drivers. Applying AlloQuant to CDK1, we demonstrate that binding of the Cyclin B1 (CCNB1) cofactor mechanically decouples a hyper-rigid inactive kinase core, allowing activating phosphorylation (pT161) to subsequently re-impose localized tension on the catalytic machinery. Conversely, the C-terminal Src kinase (CSK) faces a distinct conformational trap. While nucleotide-free monomeric CSK spontaneously samples a pre-active geometry, ATP binding excludes the active C-In conformation in all but 1 of 225 models. We show that docking partner engagement overcomes this blockade. Autophosphorylation of SRC at the activation loop (Y419) redistributes SRC conformational states without altering bulk rigidity. This redistribution is structurally coupled to the conformational state of CSK via the regulatory spine, not the catalytic machinery. Rather than mechanically deforming CSK, SRC engagement acts by conformational selection, committing roughly a quarter of CSK molecules to a fully active state. Thus, trans-allosteric kinase activation operates by defining the accessible conformational landscape of the receiver kinase. That control is exerted through mechanical remodeling in cofactor-dependent complexes and through conformational selection in transient kinase-kinase heterodimers. These findings establish AlloQuant as a general framework for quantifying how a binding partner reshapes a kinase's conformational landscape, applicable across the kinome because it assigns landmarks by profile-HMM alignment.